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ABSTRACT 

Tidal forces acting on galaxies in clusters lead to a strong dynamical evolution. In order to 
quantify the amount of evolution, I run self-consistent iV-body simulations of disk galaxies for a 
variety of models in the hierarchically forming clusters. The tidal field along the galactic orbits 
is extracted from the simulations of cluster formation in the = 1; Qq = 0.4; and f^o = 0.4, 
17a = 0.6 cosmological scenarios. For large spiral galaxies with the rotation speed of 250 km s _1 , 
tidal interactions truncate massive dark matter halos at 30 ± 6 kpc, and thicken stellar disks by 
a factor 2 to 3, increasing Toomre's parameter to Q > 2 and halting star formation. Low density 
galaxies, such as the dwarf spheroidals with the circular velocity of 20 km s _1 and the extended 
low surface brightness galaxies with the scale length of 10 — 15 kpc, are completely disrupted 
by tidal shocks. Their debris contribute to the diffuse intracluster light. The tidal effects are 
significant not only in the core but throughout the cluster and can be parametrized by the critical 
tidal density. The tidally-induced evolution results in the transformation of the infalling spirals 
into SO galaxies and in the depletion of the LSB population. In the low firj cosmological models, 
clusters form earlier and produce stronger evolution of galaxies. 



Subject headings: stellar dynamics — galaxies: clusters: general 
interactions — methods: Af-body simulations 



galaxies: evolution — galaxies: 



1. Introduction 

Hierarchical formation of clusters of galaxies in- 
volves strong local variations of the potential field. 
These tidal perturbations affect the dynamics of 
infalling disk galaxies in a number of ways. Firstly, 
they lead to the truncation of stellar and dark mat- 
ter distributions at a critical radius of the Roche 
lobe. Secondly, they increase the stellar energy 
and the amount of random motion in the disk. 
Tidal heating may be strong enough to completely 
unbind small galaxies and significantly transform 
larger ones. Also, tidal heating is likely to sta- 
bilize stellar disks against gravitational instability 
and suppress subsequent star formation. The re- 
sult of all these effects may be the morphological 
transformation of a spiral galaxy to an SO type. 

Such transformation is indeed revealed by the 
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photometric and spectroscopic studies of clusters 
of galaxies at various redshifts (Jones, Smail & 
Couch 2000, and references therein). While the 
fraction of elliptical galaxies stays roughly the 
same, the abundance of spirals decreases by a fac- 
tor of 2 to 3 from redshift 0.5 to the present. Cor- 
respondingly, the number of SO galaxies grows by 
the same amount. 

In Gnedin (2003; hereafter Paper I) I have ar- 
gued that an important role in such a transforma- 
tion is played by the tidal interactions of galaxies 
with their massive neighbors and with the global 
cluster halo. Using numerical simulations of the 
representative clusters of galaxies in three differ- 
ent cosmological scenarios, I have estimated the 
amount of tidal heating of the infalling galaxies. 
The main contribution comes from the interaction 
with the structures of order 100 kpc in size. 

Previous studies (Richstone 1976; Farouki & 
Shapiro 1981; and especially, Moore et al. 1998, 
1999; see more references in Paper I) have shown 
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that close encounters between disk galaxies can 
significantly modify their structure and leave be- 
hind the distorted spheroidals. The interaction of 
galaxies with the cluster halo has also been stud- 
ied (Merritt 1984) in the case of static spherical 
cluster models. This leads to the tidal truncation 
of galaxies. 

In this paper, I investigate the dynamical ef- 
fects of tidal heating of disk galaxies in the hier- 
archically forming clusters. I use a high-resolution 
self-consistent iV-body code to resimulate a sam- 
ple of galaxies from the lower resolution cosmolog- 
ical simulations described in Paper I. The changes 
in stellar orbits result from the effects of exter- 
nal tidal fields, computed from the second spatial 
derivative of the cluster potential along the galaxy 
trajectories. The fast variation of the tidal force 
produces tidal shocks responsible for the trunca- 
tion of galactic halos and the kinematic heating 
of stars. In agreement with the estimates of Pa- 
per I, I find stronger tidal heating in the low 
cosmological scenario. 

I use analytical arguments and discuss the ex- 
pected evolution of galaxies via tidal heating in 
§2. In §3, I describe the selection of the sample 
of galaxies from the cluster simulations. In §4, I 
present the numerical method and the tests of nu- 
merical relaxation. Then I describe the evolution 
of large spiral galaxies in §5, including the pos- 
sible effects on the gas. The disruption of dwarf 
and low surface brightness galaxies in §6 leads to 
the discussion of the origin of intraclustcr light. I 
compare these results with the previous studies in 
§7 and discuss their observational implications in 
§8- 

2. The Critical Tidal Density 

Let us start with an analytic estimate to indi- 
cate the range of systems where tidal forces are im- 
portant. The external tidal force is —(cPQ/dradr^orp = 
F a /3 rp. The amplitude of the tidal field relates di- 
rectly to the tidal stripping density, which is anal- 
ogous to the instantaneous density of the Roche 
lobe around galaxies but depends on the high- 
est tidal peaks along the galactic orbit. If a part 
of the galaxy has lower density than the critical 
value, it would be stripped or torn apart. From 
Poisson's equation, the trace of the tidal tensor 
is Ftid = F aa = 47rG j otid- A typical value of the 



tidal peaks in Paper I is F t id ~ 100 Gyr 2 , and 
the corresponding tidal density is 

^ = 1.8x10-3 ( i55 ^ =J )iif pc-». (1) 

This should be compared with the combined av- 
erage density p av = M(r)/(47rr 3 /3) of stars and 
dark matter in a galaxy. 

Consider a Milky Way-type galaxy with a large 
isothermal halo and an exponential disk (as in 
§4.1). The average halo density at the core radius 

is PhM R c) ~ °- 74o "250^k P c M © P c ~ 3 > where 
0250 is the velocity dispersion in units of 250 km 
s _1 . The average disk density at one scale length 
is pd, av {Ro) ~ °- 63 ^rf,io^kpc M © P c ~ 3 ' where 
Md,w is the disk mass in units of 10 10 M Q . Thus 
the inner regions of large spiral galaxies are stable 
against tidal disruption. 

Extended halos, falling outside the Roche lobe, 
would be truncated at a radius R t such that 
Pzv{Rt) ~ Ptid- For our models this corresponds to 
the tidal radius R t ~ 42 0250 kpc. In addition to 
the instantaneous stripping, secular tidal heating 
would enhance the mass loss. Stellar disks with 
i?o ~ 3 kpc are better protected - 99% of their 
mass is within the tidal radius. 

On other hand, low surface brightness galax- 
ies (including the dwarf spheroidals) arc in dan- 
ger of being disrupted by tidal forces. The LSB 
galaxies have lower density than normal spirals, 
< 10~ 2 Mq pc -3 , and are dominated by dark mat- 
ter at all radii. Thus their disks would hold as long 
as they are being protected by the halos. Even 
though the halo velocity dispersions are similar to 
those in large spiral galaxies, the larger sizes imply 
lower density and so stronger truncation. Stirred 
further by tidal shocks, the LSB galaxies would 
lose most of their mass, leaving behind extended 
tidal streams of stars and dark matter. These 
qualitative arguments are supported by detailed 
simulations in §5 and §6. 

2.1. Stripping of Dark Matter Halos 

The total amount of dark matter stripped from 
the galaxies can be estimated analytically in the 
spirit of Merritt (1984). Assume that all mass 
of the cluster belonged initially to the galaxies, 
distributed according to the Schechter luminosity 



2 



function 

N(L)dL = N Q ^y a e- L ' L * d ± (2) 

with a w 1. Let the mass-to-light ratio be con- 
stant and use the Faber-Jackson (1976) relation 
for the velocity dispersion, a g = cr Sj * (L/L*) 1 / 4 . 
The characteristic break of the luminosity func- 
tion, L*, corresponds to the mass scale M* w 
5 x 10 12 h~ 1 M Q . Finally, assume that the density 
profiles of the galaxies and the cluster are isother- 
mal with the three-dimensional velocity disper- 
sions u g and a c i, respectively. The characteristic 
values in the simulations of Paper I are a g = 250 
km s _1 and er c ; = \/3 x 660 km s _1 . Then the 
galactic halos will be truncated at 

Rt ~ — — R p , (3) 

where R v is the distance of closest approach to 
the cluster center (pericenter). For many galaxies 
this distance is of the order the cluster core radius, 
R p ~ 200 kpc, which gives the truncation radius 
Rt w 44 kpc. The halo mass is then reduced to 

1 cr 3 

M t tt-^-Rp^ 6.3 xlO n M . (4) 

Lr (Tel 

Tidal truncation modifies the mass function of 
galaxies such that the effective exponential cut- 
off occurs at M t instead of M*. Merritt (1984) 
estimates the fraction of mass remaining bound to 
the galaxies as 0.7 Mt/M*, or about 6% in our 
case. Tidal heating would further reduce the mass 
of the galaxies. 

For a massive cluster CL 0024+1654, Tyson, 
Kochanski & Dell'Antonio (1998) were able to in- 
fer the mass distribution using image reconstruc- 
tion of the strong gravitational lensing. The po- 
tential of the cluster appears to be very smooth, 
with only 2% of the mass bound to luminous 
galaxies. This is consistent with our estimate, as 
the velocity dispersion of CL 0024+1654 is twice 
that of our clusters. Thus even if the galaxies con- 
tained all of the cluster mass initially, they would 
retain only a small fraction of it after tidal strip- 
ping. 



3. Constructing the Sample of Galaxies 
from the Simulations of Hierarchical 
Cluster Formation 

The sample of galaxies to be re-simulated is 
chosen from the lower resolution cosmological sim- 
ulations reported in Paper I. The Particle-Mesh 
(PM) simulations use constrained initial condi- 
tions to produce representative clusters of galaxies 
for the three cosmological scenarios: Oo = 1 (clus- 
ter I), Q = 0.4 (cluster II), and Q = 0.4, Q A = 
0.6 (cluster III). All clusters have the same one- 
dimensional velocity dispersion, 660 km s~ 4 , and 
mass, 4 x 1O 14 M . The virial radii of the clus- 
ters are {1.0, 1.5, 1.4} hr 1 Mpc, respectively. The 
Hubble constant is taken to be H = 65 km s _1 
Mpc~ x (h = 0.65). The mass function of identified 
halos follows the prediction of the Press-Schechtcr 
theory for halos more massive than 3 x 10 11 M 
in Cluster I and ~ 10 12 M Q in Clusters II and III. 
Some of the lower mass halos are missing due to 
overmerging, because the PM resolution is only 60 
h^ 1 kpc. The orbits of 100 most massive galax- 
ies in each cluster are traced through the simula- 
tion starting from redshift z = 5 using a subset of 
the most bound galactic particles to calculate its 
center-of-mass motion. 

The galactic trajectories calculated in Paper I 
are unique among other cosmological simulations 
in that (i) they are real trajectories of the halos in 
the evolving cluster, and not equilibrium orbits in 
the final cluster potential; (ii) the motion of each 
identified halo is traced through the simulation at 
each time step with an unprecedented temporal 
resolution, ~ 10 7 yr; (iii) the trajectories are cal- 
culated only for 1 Hubble time, which limits the 
average number of orbital periods to 2-4; and (iv) 
the trajectories are affected by the triaxiality of 
the cluster halo, which can make initially radial 
orbits gain angular momentum and avoid the cen- 
tral regions. 

The orbital eccentricity of a galaxy is defined 
as e e (R a — Rp)/{R a + R p ), where R a and R p 
are the median apo- and pericenter distances, re- 
spectively, of the galactic orbit between z = 5 and 
z = 0. The median eccentricities for the sam- 
ple of galaxies in the three clusters are e mc d = 
{0.27, 0.53, 0.38}. They are lower than in the high- 
resolution simulation of Ghigna ct al. (1998), who 
integrated the orbits in the spherically-symmetric 
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potential approximating the cluster at the present. 
The different method of calculating orbits may be 
the cause of the difference in the results. 

In addition to the smaller eccentricities, the 
physical (not comoving) pericenter distances of 
galaxies in Cluster I are also larger. For the three 
cluster samples, i? Pimo d = {560, 160, 290} kpc. 
The pericenter distances may be relatively large 
because the clusters move throughout the simula- 
tion as they form hierarchically. Since the galaxies 
are first identified at z = 5, the cluster center of 
mass travels on average 5.2 comoving Mpc. As a 
result of cluster evolution, the galaxies often pass 
at a significantly larger distance from the center 
than they would have in a fixed cluster model. 

The rcsimulated galaxies are chosen from each 
of the three cluster samples. The galaxies ranked 
#4 by mass represent large spirals and are de- 
noted {L1,L2,L3} for the cluster models {I, II, III}, 
respectively. The orbital eccentricities of these 
galaxies, e = {0.27, 0.65, 0.38}, are typical of their 
clusters. Two additional spirals, L3a and L3b 
(ranked #13 and #73), are chosen from the Clus- 
ter III sample in order to investigate the variance 
of the results within a given cosmological model. 
In comparison with galaxy L3, they have lower ec- 
centricities (e = 0.27 and 0.09, respectively) and 
smaller pericenter distances (R p = 220 and 160 
kpc versus R p = 230 kpc for L3). 

Three low surface brightness galaxies {LSB1, 
LSB2, LSB3} arc modeled using the same orbit 
and external tidal field of galaxy L3 but with dif- 
ferent masses and sizes. Also, one galaxy from 
the bottom of each cluster sample represents the 
dwarf spheroidal galaxies, {D1,D2,D3}. Halos Dl 
and D2 are ranked #93 by mass, while halo D3 
is ranked #87. They all have low eccentricities, 
e = {0.06, 0.09, 0.24}, and orbit in the outskirts of 
their clusters. 

The initial parameters of the selected models 
are listed in Table 1. The large galaxies are fairly 
massive spirals, with the peak of circular veloc- 
ity at 250 km s _1 . The fiducial value of the halo 
mass, Mh = 8 x 10 11 M©, is not much lower than 
the true masses of L1-L3 in the cluster simula- 
tions (~2x 10 12 M Q ). The fiducial dwarf model 
had to be scaled down, because all identified ha- 
los in the cluster simulations are relatively massive 
(> 10 11 M©). This would not change the tidal his- 
tories of the dwarfs because they are still too small 



to affect the evolution of the cluster potential. 

4. Numerical Method 

The dynamical evolution of galaxies is modeled 
using the Self-Consistent Field method (Hcrnquist 
& Ostriker 1992). In this type of collisionless N- 
body code, the potential and density are expanded 
in a series of biorthogonal basis functions 

$(r) = ^A nlm <& nlm {v), (5) 
p(r) = y^A nim Pnimjr), (6) 

which satisfy individually Poisson's equation: 

V 2 $„ (m (r) = 4irGp n i m (r). (7) 

This assures that the linear combinations of the 
basis functions satisfy Poisson's equation auto- 
matically. The basis functions employed by Hern- 
quist & Ostriker (1992) consist of the ultraspher- 
ical polynomials of r and the spherical harmonics 
Yi m (6,4>). The expansion is truncated at 
10 and l max — 6, with m running as usual from — / 
to I. 

The SCF code has been modified for the current 
problem by Gnedin & Ostriker (1999) and Gnedin, 
Hernquist & Ostriker (1999). It has been opti- 
mized to run on parallel computers and offers an 
almost ideal scaling with the number of CPU. The 
simulations were done on the SGI Origin 2000 su- 
percomputer, using 10 6 particles each in the disk 
and in the halo. The simulation with 2 x 10 6 parti- 
cles takes about 20 hours on 16 R10000 processors. 

The capability to use a large number of parti- 
cles in the parallel SCF code offers a high mass 
resolution of the galactic simulations. The masses 
of stellar particles in the models of large spirals 
are 4 x 10 4 M© and in the models of dwarfs only 
100 Mq. Halo particles are more massive, how- 
ever, in proportion to the total halo mass. Also, 
the simulations gain from the capability to resolve 
temporal variations of the external tidal field on a 
scale of 10 7 yr. 

All simulations start with a galaxy running in 
isolation for one rotational period T rot at the half- 
mass radius of the disk (R\/2 ~ l-7.Ro)- The par- 
ticles reach dynamical equilibrium but do not ex- 
perience yet any appreciable heating. Then the 
external tidal field is applied for the amount of 
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time elapsed since z = 5 until z = in the three 
cosmological scenarios (T s i m in Table 3). 

The computations arc done in dimcnsionless 
units, G = Md = Ro = 1- Since the components 
of the tidal tensor F a p have dimensions of Gyr~ 2 , 
they are converted to the internal units using only 
one parameter, to = (GMd/ R^)^ 1 ^ 2 . For exam- 
ple, for the models of large spirals the conversion 
factor is t = 1.2 x 10 7 yr, or about 10% of the 
rotational period. 

4.1. Initial Galactic Models 

The initial conditions for the disk galaxy simu- 
lations were generated using a set of routines from 
Hernquist (1993). The stellar disk is a double- 
exponential, with the density profile 



Pd(R,z) = 



M d 



AttRI 



e -R/Ro e - 



\z\/z 



z 



(8) 



where Md is the disk mass, i? is the radial 
scale length, and zq is the vertical scale height. 
The dark halo is spherically symmetric and corre- 
sponds to the isothermal sphere over some radial 
interval: 



Ph{r) = 



Mi, 



2tt 3 / 2 R t (r 2 + R 2 ) 



2 IK 



(9) 



where Mh is the halo mass, R c is the core radius, 
and R t is the initial truncation radius. The nor- 
malization constant a is 



(l-^e« 2 [l-er%)]) \ 



(10) 



where q = R c /R t . Halos are effectively truncated 
at R t to reduce the computational task of follow- 
ing the orbits of very distant particles. In simu- 
lations with the external tidal field, the halos are 
always stripped to a smaller radius so that the ini- 
tial choice of R t is unimportant. For models Ll- 
L3 and D1-D3 I take Rt = 20 R c , and for mod- 
els LSB1-LSB3 R t = 6R C . In all initial models 
Rc = Ro- 

Figure 1 shows the disk and halo contributions 
to the circular velocity 



V c (r) = 



GM(r) 



1/2 



(11) 



of the initial models of large spirals. The disk dom- 
inates over the inner scale length, R , and the halo 



dominates at larger radii. The combined rotation 
curve is flat out to 10i?o- 

The figure of the halo is determined using 
the eigenvalues of the reduced tensor of inertia, 
lij = J XiXjdM. The particle coordinates are ro- 
tated into a new frame where the inertia tensor 
is diagonal. The axis ratios of the triaxial el- 
lipsoid are calculated as b/a = {I22/ hi) 1 ^ 2 and 
c/a = (I33/I11) 1 / 2 , with c < b < a. The halos are 
initially spherical with an isotropic velocity distri- 
bution. 

The initial orientation of the galactic disk with 
respect to the coordinate frame of the cluster simu- 
lation from Paper I can be set arbitrarily. I assume 
that the axes of the large galaxies are parallel to 
the cluster frame and the axes of the dwarf galax- 
ies are tilted by 45°. The tidal tensor has compa- 
rable components in all directions, and therefore 
the results should not depend on the initial disk 
orientation. 

4.2. Tests in Isolation 

Before studying the effects of tidal heating, it 
is very important to assess the extent to which 
the simulations are affected by numerical artifacts. 
The main is numerical relaxation, which results 
in erroneous heating of stellar orbits. This re- 
laxation arises from Poisson fluctuations of the 
coarse-grained potential represented by a finite 
number of particles, N ~ 10 6 . This is much 
smaller than the number of stars in real galax- 
ies, and therefore the relaxation rate is artificially 
high (e.g. Weinberg 1996). 

I have run five test models in isolation in order 
to check the dependence of numerical heating on 
the number of particles and the time step (Tables 
2a and 2b). The first model, L-isol, is the main 
model for large spiral galaxies L1-L3 and the next 
three are its modifications with reduced number 
of particles or increased time step. The last iso- 
lated model, LSB-isol3, is that of a large low sur- 
face brightness galaxy LSB3. The tests run for 
12 Gyr, the longest simulation time of the three 
cosmological scenarios. 

Due to the diffusive nature of the relaxation 
process, both the first and second order changes 
of stellar energy are comparable, (AE/E) <~ 
((AE/E) 2 y However, the former may change 
sign and vary randomly, whereas the latter rises 
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monotonically in time and offers a better mea- 
sure of relaxation. The average energy change of 
disk particles at the end of run L-isol is small, 
((AE/E) 2 ) w 0.5%. Table 2a shows that the 
radial mass distribution of the disk remains ex- 
ponential and the scale length, R , changes only 
by 4%. The halo also remains spherical, with the 
axis ratios b/a and c/a close to unity. The only 
noticeable difference is the vertical heating of the 
disk. The average scale height, zq, increases by a 
factor 1.7 and the vertical velocity dispersion a 2 z 
rises similarly. As Figure 6 shows (lower dots), the 
thickening is stronger near the center where the 
density of particles is higher and the relaxation is 
expected to be faster. 

Another statistic of disk heating is the change 
of the linear vertical energy. In equilibrium it can 
be defined using the epicycle theory as 



E z 



2 + 2U 2 Jo 



(12) 



where the second derivative of the potential is eval- 
uated in the plane of the disk (for a thin disk it is 
2ttGY^(R)/z , where £(i?) is the surface density). 
Equation (12) is strictly an integral of motion only 
near the midplane, where the disk dominates grav- 
ity. In the isolated runs the vertical energy E z 
changes significantly more than the total energy 
E, by the order itself (Table 2b). But this is con- 
sistent with the changes of the scale height and 
the vertical velocity dispersion. 

With half the number of particles in model 
L-isol2, the second order energy change is 0.8%. 
The vertical energy and the scale height are pro- 
portionately larger, while the radial scale length 
and the figure of the halo are unchanged. If only 
the number of disk particles is reduced (model L- 
isol3), the heating is slightly less but similar to the 
previous case. Thus only the vertical structure of 
the disk is sensitive to numerical relaxation, when 
the number of particles in simulation N <~ 10 6 . 

The time step in the main simulations is At = 
0.1 to in code units, or about 1% of the half-mass 
rotation period, T ro t = 27ri? 1 / 2 /V r c . In a thin disk 
the stellar orbits should be calculated accurately 
even in the center, where the rotation period de- 
creases by 60% as the disk density increases by a 
factor of 5. However, once the relaxation raises 
velocities artificially high, the stellar motion devi- 
ates more and more significantly from the original 



orbits. Doubling the time step in model L-isoll 
leads to larger energy changes than in the other 
isolated models. This effect is coupled to numer- 
ical relaxation and, therefore, both are more im- 
portant in the inner parts of the galaxy. This is 
not a problem for the present study because, by 
construction, external tidal forces vanish in the 
center. Tidal effects are more significant in the 
outer regions of the galaxy where both numerical 
artifacts are less important. 

While the artificial thickening of the disk 
is certainly undesirable, it can be corrected 
for in a statistical way. The relaxation rate, 
d ((AE/E) 2 ') /dt, is constant in time and the scale 
height zo,num increases linearly (cf. upper panel of 
Figure 10), as expected for pure Poisson fluctua- 
tions. Thus at any time it can be subtracted from 
the real simulation to accentuate the effects of 
tidal heating. Assuming that the numerical heat- 
ing is statistically independent from the external 
tidal heating, the two effects add in quadrature 
and the modified scale height can be defined as 

z~o = (zq — z 2 num + 2q i ) 1 / 2 , (13) 

where zqa is the initial scale height. The least 
squares fit to the growth rate is 

dlnz , nU m/dt = (6.0 ±0.05) x 10~ 2 Gyr" 1 . (14) 

In the following section we derive the rate of nu- 
merical relaxation in a thin disk to confirm that 
the growth of zo,num(£) is consistent with the self- 
interaction of disk particles and with the heating 
by more massive halo particles. 
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4.3. Numerical Relaxation in a Thin Disk 



Spitzer & Hart (1971) have derived the median relaxation time for a smooth spherical system of mass M 
and half-mass radius Rh'. 



trh = 0.U-^-(^-) . (1.1) 



In A \GM , 

It describes very well the rate of heating in simulations, even including the resonant interactions of stars (see 
the discussion in Gncdin & Zhao 2002, based on the results of Weinberg 1993). The relaxation is dominated 
by distant encounters with large impact parameters b max , parametrized by the Coulomb logarithm, In A = 
hifrmax/frmin- The system relaxes on many (~ 0.1N) dynamical times and for large enough TV can be 
considered collisionless. 

However, Rybicki (1971) has emphasized that the relaxation time is much shorter in the infinitely thin 
disks, of the order an orbital period (see also Binney & Tremaine 1987, problem 8-6). This is due to the 
close encounters dominating in 2D geometry, instead of the distant encounters in 3D geometry. In a disk of 
finite thickness this effect largely disappears, but there is another effect that leads to faster relaxation. The 
collective response of stars to a perturbation is greatly enhanced when the stellar motions are highly ordered 
in a kinematically cold disk (Sellwood 1987). Since the stellar orbits are still largely confined to a plane, the 
most noticeable effect of disk relaxation is the heating in vertical direction. 

4- 3.1. Heating by Disk Particles 

For the general case of a thin disk with an initial scale height zq^ <C Rq and vertical velocity dispersion 
&z,i *C V c , the relaxation rate can be estimated as follows. The velocity kick in the vertical direction from a 
single encounter with a star m* at an impact parameter b and relative velocity v is Av z ss 2Gm*/vb. The 
encounter velocity is of the order the radial and azimuthal dispersion in the plane. Most of stellar encounters 
take place at a median height above the plane, where the density of the double-exponential disk (eq. [8]) is 
half of that in the plane, p d (R,zi/ 2 ) = 1/2 x Y,{R)/2zq. The integrated effect of encounters at all impact 
parameters is then 

A min V v b ) 42om* vzo(t) 

where 6 max ~ 2zq ~ 0.6 kpc and 6 m i n ~ 0.1 kpc (set by the resolution scale). In the simplest case, as a result 
of heating the surface density S(i?) remains the same but the scale height zq increases in time. The rate of 
increase can be calculated from the relation a\ = ttGT,(R)z , which should hold as long as the disk is still 
relatively thin: 

dz _ 1 da\ _ 2Gm»lnA d 

~dt ~ 7rGS(i?) ~dt ~ ^ ' ( ' 

The solution is 



zo /, 2* N 1/2 



where the vertical relaxation time is 



%0A V ^rel,: 



N d vzl 



(18) 



^^cwr 11 ^^ (19) 

for v = 80 km s" 1 , M d = N d m* = 4 x 10 10 M Q , z ^ = 0.3 kpc, lnA d = 1.8, and N dfi = N d /10 6 . The 
corresponding rise of the vertical velocity dispersion is 

^o\ 1/2 A , 2t ^ 1 / 4 



. • i 1 + 7 — I • iW) 
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This matches qualitatively the heating of the disk by massive molecular clouds (cf. Lacey 1984). The initial 
growth is linear and the maximum rate is d\nzo/dt = t~ e ] z w 9 x 10~ 2 Gyr -1 , which is in reasonable 
agreement with the observed rate of growth in the isolated simulation (eq. [14]). 

If the stars were distributed spherically rather than in the disk, the relaxation time would be significantly 
longer (eq. [15]), t r h ~ 10 12 yr. Alternatively, if the stars were confined to an infinitely thin disk, the heating 
rate (in the plane) would be 

da 2 /c >^-™ N 2 



_R — 
dt 



rf^ya M „?«. (2D 

Amin V Vb J m * Vb min 



This rate (clearly dominated by close encounters) is faster than the finite-thickness rate (dominated by 
distant encounters, eq. [16]) by a factor 2.1 z (t)/z ,i- So as expected, the finite-thickness disk provides an 
intermediate case between the 2D and the 3D geometry. 

4-3.2. Heating by Halo Particles 

When the disk is surrounded by a dark halo, the heating by massive halo particles can be more important. 
The heating rate by the particles with mass m,h and density ph is 

da 2 /-6„ax / 2Gm h \ 2 p h 8irG 2 m h p h In A h 

' -v2-Kbdb= , (22) 



J brain V V t> J m h V 

where the encounter velocity v is dominated by the halo dispersion, a h . The ratio of the heating rates due 
to the halo and disk particles is 

(do- 2 z /dt) h = m h p h InAfefcncd (2 „s 
(da 2 /dt) d m* p d lnA d w onCjh ' 

The typical parameters in our models at i?o are to^/to* = 20, v nc,h/wcnc,d ~ 2, and Pd/ph ~ 5 — 10, so 
the two heating rates are comparable. However, the growth of the scale height is different because now the 
right-hand side of equation (22) is independent of a z and 



%0,i \0~z,i J ^rel,h 

with 



(24) 



Uchh= 8n\nA h G 2 M h p h (25) 

This interaction with the halo particles is equivalent to the heating by massive black holes (cf. Lacey & 
Ostriker 1985). For M h = 8 x 10 11 M Q , v w 160 km s^ 1 , a z . t w 80 km s" 1 , In A h w 2, and p h {R c ) w 0.08 M 
pc~ 3 , the relaxation time is £ re i,h ~ 17^,6 Gyr. This is in excellent agreement with the growth rate in the 
isolated simulation (eq. [14]). Thus, both disk and halo particles contribute to the disk thickening, although 
in the long run the heating by halo particles would dominate. 

Having a certain understanding of the numerical relaxation in the disk, when it sets in (eqs. [19] and [25]) 
and how to account for it (eq. [13]), we move to the simulation results. 



8 



5. The Transformation of Large Spiral 
Galaxies 

We start with the results for large disk galax- 
ies. Figure 2 shows the distribution of particles 
at the beginning and the end of the simulation of 
galaxy LI. The first important result is that the 
galaxy survives tidal heating. It still looks like a 
disk galaxy of roughly the same size. Large, co- 
herent spiral arms form in the outer regions. The 
disk shrinks in the radial direction and flares at 
the fuzzy edges. But the main difference from the 
initial configuration is a substantial thickening of 
the disk. Viewed at a certain surface brightness 
limit, this thick and round disk galaxy may be 
classified as an SO. 

The dark matter distribution is relatively sim- 
ple. The galaxy has been stripped of about 40% 
of its halo (Table 3) and has shrunk to the limits 
of the stellar disk. Inside the truncation radius, 
the halo has been little affected (Figure 3). The 
figure of the halo becomes triaxial but only mildly, 
with the smallest axis ratio c/a = 0.88. Thus, the 
amount of the remaining halo mass is indicated 
simply by the truncation radius, Rt- For galaxy 
LI it changes from 60 kpc to 37 kpc. Figure 4 
shows that the rotation curve is still flat up to six 
disk scale lengths. The circular velocity, measured 
at 5Ro, is decreased only by 2%. 

The disk of the galaxy is perturbed more sig- 
nificantly. Figure 5 shows that the surface density 
decreases in the center but increases between 3i?o 
and 4i?o, relative to the initial profile. The surface 
density exhibits noticeable wiggles superimposed 
on the otherwise exponential disk. 

In addition to truncating the halo, a major ef- 
fect of tidal heating is the thickening of the disk. 
By construction, the external tidal force vanishes 
at the center and increases linearly with radius. 
Figure 6 demonstrates that the disk flares in the 
outer parts. Numerical relaxation affects the disk 
structure only within the inner two scale lengths, 
where tidal heating is negligible. At R > 2Ro, the 
external tidal effects are clearly dominant. 

A direct evidence of tidal heating is the dra- 
matic increase of all three velocity components in 
the outer regions of the disk (Figure 7) . The ve- 
locity dispersion increases at the expense of the 
rotation speed, which starts to decline at around 
four scale lengths from the center (while the circu- 



lar velocity is still flat; cf. Figure 4). In contrast, 
the velocity dispersion of the halo decreases in the 
outer parts due to the significant mass loss. 

The analysis of the other large spirals shows a 
lot of similarity but with quantitative differences. 
Since the tidal field is stronger in the low f2 clus- 
ters (see §5.3), galaxies L2 and L3 are subject to 
stronger perturbations than galaxy LI. Galaxy 
L2 suffers the biggest mass loss, about 64%, while 
galaxy L3 provides the intermediate case with 
about 50%. Their truncation radii fall below 30 
kpc. The disks of both galaxies L2 and L3 are 
heated more strongly than in galaxy LI. 

Figure 8 shows that the mass loss from the 
three halos is not steady. There are several in- 
stants of large losses associated with the maxima 
of the tidal force. Examining the particle distri- 
bution at successive epochs reveals how the halo 
is being stripped. After a strong tidal shock, a 
cloud of dark matter particles disperses across the 
tidal boundary and slowly drifts away from the 
galaxy in two opposite directions along the same 
line. That line corresponds to the direction of 
the strongest component of the tidal tensor at the 
time. The unbound particles are later stretched 
along that direction further and further away from 
the galaxy. 

As Figure 9 shows, the disk scale lengths of all 
three galaxies remain almost constant. In con- 
trast, the vertical scale heights rise monotonically 
with time (Figure 10). While a part of this effect 
is due to the numerical relaxation (eq. [13]), the 
corrected scale height zq measures the thickening 
only due to the tidal heating. Galaxy L2 leads 
in the rate of expansion, in agreement with other 
indicators of evolution, the mass loss and the ve- 
locity dispersion increase. At the end of the run, 
its scale height is more than doubled. The other 
two galaxies experience more modest evolution. In 
particular, the disk of galaxy L3 had expanded to 
its maximum height in 7 Gyr and then remained 
effectively unchanged for the same amount of time 
until the present. 

5.1. Dependence on Time Step and Num- 
ber of Particles 

In order to quantify the dependence of the re- 
sults on numerical details, such as the time step 
and the number of particles, I have run three ad- 
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ditional models of galaxy L3 (Table 4). They are 
simulated using the same parameters as the iso- 
lated models {L-isoll, L-isol2, L-isol3} but with 
the external tidal field of galaxy L3. The corre- 
sponding isolated runs are used to correct the scale 
height z . 

In model 1, the time step is doubled while 
the number of particles is kept the same, rela- 
tive to the main model of L3. As discussed in 
§4.2, the time step is chosen small enough to rep- 
resent accurately stellar orbits even in the center; 
At = 0.01 T rot w 10 6 yr. Indeed, with the dou- 
bled time step model 1 reproduces most of the 
features of model L3. The main changes are the 
13% stronger thickening of the disk and a slightly 
longer radial scale length. 

In model 2, half the number of particles is used 
in both the disk and the halo with the original 
time step. The results are very close to those of 
model L3. The changes of the scale length and the 
corrected scale height are not noticeable at all. 

In model 3, where only the number of disk parti- 
cles is reduced by half, the scale height grows even 
less than in model L3. This indicates that most of 
the numerical heating of disk particles is produced 
by the interaction with massive halo particles, as 
expected. In all three test runs the halo structure 
and mass loss are virtually identical. 

The conclusion from these tests is that only the 
vertical structure of the disk is sensitive to nu- 
merical effects. Variations of the time step and 
the number of particles affect the corrected scale 
height z~o by 10% to 30%. Thus the deviation of 
our main results from the case of infinite resolution 
is expected to be of the same order. 

5.2. Dependence on Force Resolution 

The necessary truncation of the spherical har- 
monics series in equation (5) imposes a limit on 
the vertical force resolution in a thin disk. An ac- 
curate calculation of the force with z /R = 0.1 
may require the expansion order ? max ~ 100. 
Our simulations with Z max = 6 thus have an ef- 
fectively smoothed vertical force within the disk, 
and the comparison with the forces calculated by 
direct summation shows significant deviations at 
\z\ < 0.4 Rq. From the SCF simulations alone it 
is difficult to assess how this could affect the re- 
sponse of the disk to external perturbations. In or- 



der to test sensitivity of the results to force resolu- 
tion, I use the publicly-available tree code gadget 
(Springcl, Yoshida & White 2001) and resimulate 
large galaxy models LI, L2, and L3. 

The tree algorithm with an accuracy parameter 
6 = 0.8 calculates the vertical forces in our thin 
disks within 1% of the direct summation result. 
The algorithm is, however, much more compu- 
tationally expensive and can be realistically used 
only with a subset of particles. With Nd = 10 5 , 
Nh = 4 x 10 5 , one parallel gadget simulation 
takes 150 hours on 8 processors of the beowulf PC 
cluster at STScI. I use a 'type 0' time step criterion 
based on local particle acceleration a and soften- 
ing parameter e: At oc \fe]~a, and limit the time 
step not to exceed 0.1 in code units (w 10 6 yr, the 
same step as in the SCF runs). The average time 
step was 0.04— 0.08 in code units. In order to ad- 
equately resolve the vertical structure of the disk 
and to minimize relaxation effects, I take the soft- 
ening parameter a = 0.1z for disk particles and 
€h = 0.1i?o for halo particles, using the definition 
of Springel, Yoshida & White (2001). 

Since the numerical heating of the disk is largely 
due to the massive halo particles, I run two sets 
of models increasing Nh from 10 5 to 4 x 10 5 . The 
number of disk particles is fixed at Nd = 10 5 (these 
numbers are limited by a realistic run time). In 
the control run in isolation with Nh — 10 5 the 
disk size shows a large increase (zo/zo,i = 3.0) but 
so does a corresponding SCF run with the same 
number of particles (zq/zqj = 3.1). In the control 
run with Nh = 4 x 10 the numerical heating is 
already reduced (zo/zo,i = 1.9) almost to the level 
of the main SCF runs. In both isolated gadget 
simulations the radial scale length changes by less 
than 1%. This agrees with the discussion in §4.3 
that numerical heating depends primarily on the 
number and masses of particles. 

Table 5 shows the results of gadget simula- 
tions for models L1-L3. Overall, the disk heating 
is somewhat stronger than in the SCF runs: both 
z and R increase by roughly 10%. The vertical 
heating is much stronger in the runs with fewer 
halo particles - this is clear evidence that numeri- 
cal relaxation can affect the evolution of thin disks 
even with adequate force resolution. 

The thickening of the disk is less monotonic in 
the GADGET runs. Sudden jumps of zq, for exam- 
ple at 5 Gyr and 9 Gyr for model L2, coincide with 
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larger than usual deviations of the disk center-of- 
mass from the halo center-of-mass. The two can 
differ by ARc m < 800 pc after a strong pertur- 
bation, and in such an unrelaxed state the disk 
appears thicker. Note that this deviation is still 
within 0.3 Rq. In models LI and L3, and in all 
SCF runs, AR Co m < 100 - 200 pc. 

The disk appears symmetric in the vertical di- 
rection and shows no sign of warping. The ori- 
entations of the inner and outer parts of the disk, 
divided at 2.5 R , differ by less than 5° (8° for SCF 
runs). 

A noticeably difference in the gadget runs is 
the radial expansion of the disks. In the SCF mod- 
els, Ro was essentially unchanged. This might be 
a result of the restoring forces being effectively 
"anisotropic" within the disk: while the radial ex- 
pansion is sufficiently accurate, the vertical expan- 
sion is smoothed. Fortunately, the difference with 
the gadget results is only 10 - 15%. 

As a result of the increase of Ro, the circular 
velocity drops correspondingly by about 10% in 
all GADGET models. However, the halo structure 
in the gadget runs almost exactly mirrors that 
in the SCF runs. The only deviation of the halo 
mass is in model L3 (by 3%), while the shapes of 
the halos are similar. I conclude that apart from 
Ro, the only quantity susceptible to the errors in 
force resolution is zq, and while it increases by a 
factor 1.6 — 2.4 in the SCF runs the true value is 
probably in the range 2 — 3. 

5.3. Dependence on Tidal History 

In cluster simulations presented in Paper I, the 
amount of tidal heating has been estimated us- 
ing the tidal heating parameter, J t id- It is sim- 
ply the square of the integrated velocity change 
summed over all peaks of the tidal force F a p 
{a,P= {x,y,z}): 

n a. ^ J / n \ dyn / 

(26) 

It is also corrected for the conservation of adia- 
batic invariants of stellar orbits when the external 
perturbations are long, using the parametrization 
of Gnedin & Ostriker (1999). Here r„ is the ef- 
fective duration of peak n for each value of a and 
P, and tdyn is the half-mass dynamical time of the 



galaxy (<~ T rot ). The average energy change due 
to tidal heating for the stars at radius r is 

(AE) = ±I tid (r 2 ). (27) 

This semi-analytical estimate provides a good de- 
scription of the energy changes of particles in the 
simulation and can be used to compare tidal ef- 
fects for different galaxies. 

The galaxies in each cluster sample show an 
exponential distribution of the tidal parameter, 
N(I t id) = N cxp(— Itid/Io)- The characteristic 
parameter Iq is the largest in Cluster II (5 x 
10 3 Gyr -2 ), intermediate in Cluster III (2 x 10 3 
Gyr~ 2 ), and the smallest in Cluster I (8 x 10 2 
Gyr~ 2 ; see Paper I). The simulations of galax- 
ies L1-L3 show noticeable differences in the evo- 
lution of large spiral galaxies in the three cosmo- 
logical models. In agreement with the estimates of 
the tidal heating parameter, galaxy L2 experiences 
stronger evolution than galaxies LI and L3. 

However, these differences may not serve as a 
clear discriminant among the cosmological mod- 
els. The distribution of tidal parameters within 
the same cluster leads to a comparable variance 
of the results. The additional simulations of two 
galaxies from Cluster III prove this point. 

The two computed models, L3a and L3b, use 
the same initial conditions as galaxy L3 but dif- 
ferent tidal histories. These galaxies, correspond- 
ing to the Cluster III halos #13 and #73, have 
some of the largest heating parameters, 7 t id = 
3.8 x 10 3 Gyr- 2 and I tid = 6.4 x 10 3 Gyr" 2 , re- 
spectively. For comparison, the reference galaxy 
L3 has Jtid = 4.2 x 10 3 Gyr -2 . These models rep- 
resent the strongest cases of tidal evolution. 

As expected, galaxy L3b loses more mass and 
develops a thicker disk than model L3, while for 
L3a the trend is the opposite (Table 3). The final 
state of galaxy L3b is similar to that of model 
L2 from Cluster II. On the other hand, galaxy 
L3a is closer to model LI from Cluster I. Thus 
the distribution of tidal histories of galaxies within 
the same cluster leads to a similar variation of the 
dynamical evolution. 

The adiabatic correction factor in equation (26) 
is very important. It accounts for the temporal 
structure of the tidal field and suppresses the per- 
turbations with long timescales to which galax- 
ies respond adiabatically. Without that factor, 
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galaxy L3a would have had a higher value of itid 
than galaxy L3b. The difference is due to the tidal 
force being more "impulsive" for galaxy L3b, while 
being more "adiabatic" for L3a. 

5.4. Do spirals become SO galaxies? 

The simulations of large spiral galaxies show 
that the kinematic heating and vertical expansion 
of the disks lead to a significant morphological 
transformation. An initially normal spiral may be 
identified in the end as a lenticular galaxy, an SO. 
The following section puts this claim to a thorough 
observational check. 

The final snapshot of galaxy LI (Figure 2) looks 
similar to the deep exposures of the SO galaxy 
NGC 4762 (Sandage 1961; Burstein 1979). The 
disk is thick and ends with fuzzy warps. A mod- 
est central mass concentration is seen in both the 
simulation and the real galaxy. 

The most detailed data are available for nearby 
SO galaxies in the field (Kent 1985). The HST 
images of lenticulars in distant clusters (Ellis et 
al. 1997; van Dokkum et al. 1998) usually provide 
only integrated light and colors. Typically, the 
surface brightness profiles of SOs are well fitted 
by a combination of the exponential disk and the 
r 1 / 4 bulge. The simulations reported here do not 
include a bulge component initially. In fact, one of 
the goals of this study was to investigate whether 
a bulge forms as a result of tidal heating. Two 
competing effects play a role here: stars on aver- 
age gain energy from tidal perturbations and move 
generally outward; on the other hand, the tidally- 
induced energy dispersion contributes to two-body 
relaxation and accretion onto the center. Figure 
5 shows that no appreciable bulge forms from the 
disk particles, asserting that the relaxation time 
is still much longer than the age of the galaxy. 
However, the interstellar gas may be driven in- 
ward by tidal effects and ignite a burst of star for- 
mation (Hernquist & Mihos 1995; see also §5.5). 
This would assist in the formation of a prominent 
bulge. 

SO galaxies typically have the same disk scale 
lengths as the spirals of similar luminosity (Kent 
1985). In the present simulations the scale length 
i?o is virtually unchanged. Also, the axis ratios 
of the simulated disks, 0.15 to 0.25, are consis- 
tent with the average observed value of 0.2-0.25 



for the "thick disks" (Binney & Merrifield 1998). 
Burstein (1979) attempted to decompose the verti- 
cal brightness profile of late type edge-on galaxies 
into a combination of the thin and thick disks. The 
latter dominate in the outer regions of the galax- 
ies and are much more prominent in SOs than in 
spirals. Thick disks are not as centrally peaked as 
thin disks and have a shallow radial gradient: they 
look "diffuse" and "boxy". All this is consistent 
with having a tidal origin; flaring of the disk was 
expected and has been detected in the simulations. 

Kinematic information is available only for the 
local galaxies. In a cold exponential disk, vertical 
velocity dispersion should fall with radius as the 
surface density: a^{R) = 2ttGzo^o exp(— R/Rq). 
In a sample of ten SO galaxies (Seifert & Scorza 
1996) the velocity dispersion flattens at 50 — 100 
km s _1 and even rises in the outer regions. Sim- 
ilarly, a large study of the early type galaxies by 
Simien & Prugniel (1998) shows a systematic in- 
crease of the total velocity dispersion at large radii, 
often exceeding 100 km s _1 . Present simulations 
certainly confirm this trend. 

Finally, note that bars are commonly found in 
SO galaxies. A single tidal interaction may excite 
bar instability in an axisymmctric galaxy, but a 
succession of tidal shocks and the resulting kine- 
matic heating would isotropize stellar orbits and 
destroy the bar. The actual survival of bars, then, 
depends on the amount of tidal heating in the in- 
ner galaxy (see Sellwood 1999). Present simula- 
tions, which include a stabilizing dark matter halo, 
do not show any detectable bar. 

5.5. Gas Dynamics and Star Formation 

The evolution of spiral galaxies depends sensi- 
tively on the amount and distribution of the in- 
terstellar gas. Even though the collisionless sim- 
ulations do not address gas dynamics directly, it 
is possible to infer qualitative features of the gas 
response to tidal heating. 

The external tidal force is not steady in time, 
it comes in a series of peaks (cf. Figures 14-16 in 
Paper I) . The timescale of each tidal perturbation 
is of the order 10 s yr or longer. While the stars 
would continue to gain energy without dissipation, 
the gas should cool to its equilibrium temperature 
fairly quickly. The cooling time of the ionized gas 
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gas, (5, required to make the disk unstable. Writ- 
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where A(T) ~ 10 -23 erg cm 3 s _1 is the cooling 
function at the temperatures T ~ 10 6 K, and n 
and n e are the baryon and electron number densi- 
ties, respectively. As a case of extremely fast tidal 
heating, consider that the gas with nT ~ 10 4 K 
cm~ 3 heats up to 8 x 10 5 K, corresponding to the 
velocity dispersion of 100 km s -1 . Even in this 
case, it would cool back to 10 4 K in 4 x 10 7 yr. 
In reality the gas would likely heat less efficiently 
and therefore cool faster. 

However, the thick stellar disk is now stable 
against gravitational perturbations and is unlikely 
to form new stars. Figure 11 shows Toomre's 
Q parameter (Toomre 1964) for the final states 
of galaxies L1-L3. Even the minimum value, 
Qmin ~ 2, is large enough to suppress any axisym- 
metric disk instability. Under such conditions spi- 
ral structure and gaseous shocks would not form 
in the inner regions of the galaxies. The spiral 
arms visible in Figure 2 are prominent only in the 
outer regions, where the diffuse gas is likely to be 
ionized or stripped by the intraclustcr medium. 

The stability of stellar disks is modified in the 
presence of cold gas. The low velocity dispersion 
of the gas amplifies the growth of perturbations 
and lowers the effective value of Q compared to 
the pure stellar case. Jog & Solomon (1984) de- 
rive a stability criterion for a two-fluid system of 
the stars and gas. From their equation (22) it fol- 
lows that the system is unstable if the modified 
parameter Q c s is less than unity: 



Qeff = 



3.36 
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<1, (29) 



where [3 — £ gas /E* is the ratio of the gas to stellar 
surface densities, c = vth/c* is the ratio of the gas 
thermal velocity to the stellar velocity dispersion, 
and x = k/(ct* k) is the ratio of the perturbation 
wavelength, 2-K/k, to the effective scale of stellar 
motion. Here 
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is the epicyclic frequency, which includes the con- 
tributions from both stars and dark matter. Equa- 
tion (29) allows one to determine the amount of 



ing the velocity dispersions as u t h = H T £ 
km s _1 and cr» = 
c = 0.11 T 4 1/2 a-} 00 
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100 cvioo km s 1 , I find 
The most unstable mode is 
usually a half of the maximum allowed wavelength 
(Binney & Tremaine 1987), which in this case is 
Xmax ~ 1.87(1 + (i)/Q. With these parameters, 
the equation Qcff = 1 can be solved numerically 
to yield (3 ~ 0.52. It shows that the disk is unsta- 
ble only when the mass of cold gas is at least 50% 
of the stellar mass, a condition rarely satisfied in 
any spiral or SO galaxy even in the field. 

Current theories of star formation in normal 
spiral galaxies often rely on a continuous infall of 
cold gas (e.g., Edvardsson et al. 1993; Dalcanton 
et al. 1997b). This infall supplies the material for 
new stars and also helps to sustain a marginal dy- 
namical instability of the disk. After galaxies enter 
the cluster, external tidal forces halt the late infall 
and, with it, new star formation. 

5.6. Ram Pressure Stripping 

The amount of gas remaining in the galaxies 
can also be reduced by a variety of processes. The 
intracluster media (ICM) can sweep the diffuse gas 
via ram pressure (Gunn & Gott 1972). Assuming 
the ICM density in the cluster core nj cm = 10~ 3 
cm -3 (Jones & Forman 1984), the surface density 
of H I gas E gas = 4 M Q pc~ 2 (Binney & Mer- 
rifield 1998), and the orbital velocity of 800 km 
s , large spirals with exponential disks would be 
stripped of atomic gas down to R = 5 kpc. Recent 
numerical simulations of ram pressure stripping in 
rich clusters (Abadi, Moore & Bower 1999) con- 
firm that the H I disks are truncated at ~ 4 kpc. 
Through viscous coupling even the inner gas may 
be affected (Quilis, Moore & Bower 2000). 

The efficiency of ram pressure stripping de- 
pends also on the hierarchical formation of the 
cluster. When individual groups of galaxies merge 
into the cluster, a bow shock may form slightly 
ahead of the group's motion. It happens when the 
ram pressure on the infalling gas increases faster 
than a sound wave can propagate across the sub- 
structure (Fabian & Daincs 1991). Numerical sim- 
ulations (Rocttiger, Loken & Burns 1997) confirm 
that bow shocks do form around the substructures 
falling in at supersonic speeds. These shocks may 
delay the heating of gas to the virial temperature. 
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X-ray observations of the central region of Coma 
cluster by ASCA (Hank Donnelly et al. 1999) find 
a cold and a hot spots in the temperature map. 
A natural explanation would be a bow shock as- 
sociated with the group around the bright galaxy 
NGC 4869. 

While ram pressure stripping does not affect 
directly dense molecular clouds, the interpenetrat- 
ing galactic encounters may remove the molecular 
component (Kritsuk 1983). Direct collisions be- 
tween molecular clouds would be rare but the cou- 
pling via magnetic fields would ensure dynamical 
interaction. In high-speed collisions the gas would 
shock and decelerate, whereas the stars and dark 
matter fly away. 

In a survey of 17 bright spirals in the Virgo 
cluster, Cayatte et al. (1994) find two groups of 
galaxies with different H I properties. One group 
has a typical H I surface density (~ 4 M Q pc~ 2 ) in 
the inner regions, but a steeply falling one beyond 
a half of the optical radius. These galaxies could 
be affected by ram pressure stripping. The other 
group has lower central surface densities but the 
gas extends much further, often beyond the op- 
tical radius. In other words, the latter galaxies 
are modestly deficient of H I at all radii, indepen- 
dently of the local gravitational force. Other evap- 
oration mechanisms may play a role here, such 
as the turbulent mixing or heat conduction (see 
Sarazin 1988 for details). 

Also, in a study of three rich Abell clusters 
and the Virgo cluster, Valluri & Jog (1991) find 
that the H I deficiency is stronger in larger spi- 
rals. This statistical trend is opposite to that ex- 
pected for ram pressure stripping, since more mas- 
sive galaxies should hold on to their gaseous com- 
ponent better. Having considered several possible 
mechanisms, Valluri & Jog (1991) conclude that 
the H I deficiency can only be explained by tidal 
interactions. 

Finally, note that numerical simulations of 
galaxy interactions (Hernquist & Mihos 1995; 
Barnes & Hernquist 1996) often show a dramatic 
loss of angular momentum of the gas, drawing it to 
the center and igniting a powerful starburst. This 
is corroborated by the observation that galaxies 
with distorted profiles are more likely to be de- 
tected in Ha than the undistorted ones (Kcnnicutt 
1998). Moss & Whittle (1993) find a 50% higher 
detection rate of Sa galaxies in the eight nearby 



Abell clusters than in the field. These starbursts 
may also blow away the remaining gas from small 
galaxies (Dekel & Silk 1986), leaving behind only 
compact stellar spheroids. 

6. The Disruption of Low Density Galax- 
ies 

6.1. Dwarf Spheroidal Galaxies 

Dwarf galaxies are abundant in the Universe 
and dominate the number counts both in the field 
and in clusters (Ferguson & Binggcli 1994; Gal- 
lagher & Wyse 1994). While in the field they 
are usually found around large luminous galax- 
ies, their distribution in clusters is dramatically 
different. There are many fewer dwarfs in the cen- 
ter than on the periphery (Phillipps et al. 1998) 
and there is a strong anti-correlation of the ratio of 
dwarfs to giants with the local projected density of 
galaxies (Driver, Couch & Phillipps 1998). Also, 
less massive clusters are likely to have more dwarf 
galaxies than large virialized clusters. Observa- 
tional evidence suggests that dwarfs are subject 
to strong selection effects in clusters, consistent 
with tidal disruption. 

Dwarf galaxies are usually divided into two dis- 
tinct classes, the high density dE and the low sur- 
face brightness dSph (Kormendy 1985). It is the 
dark matter-dominated dwarf spheroidals that are 
likely to be disrupted. For the purpose of this 
simulation I assume that initially they were disk 
galaxies embedded in dark matter halos. The sur- 
face brightness profiles of many dSph are fairly 
good exponentials, which supports their disk ori- 
gin. The typical scale length is about 1 kpc (Fer- 
guson & Binggeli 1994), similar to dwarf irreg- 
ulars like the LMC (de Vaucouleurs & Freeman 
1972). The velocity dispersion of dSph is very 
low, typically 10 km s _1 and no more than 20 
km s _1 (Kormendy 1985). I assume the latter ex- 
treme value to investigate whether even the dens- 
est dwarf spheroidals could survive tidal heating. 

I run three dwarf models {Dl, D2, D3} with 
the tidal histories from different cluster simula- 
tions (see §3). Table 1 shows that the stellar disks 
are small and are completely dominated by dark 
matter. Still, the density of the halo cores is low, 
Ph,av ~ 5xlO~ 3 M Q pc -3 , dangerously close to the 
critical tidal density (eq. [1]). The peaks of the 
tidal force can reach as high as 200 Gyr -2 , which 
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corresponds to the tidal density of 3.6 x 10 3 M 
pc~ 3 . 

None of the dwarfs survives until the end of 
the simulation. Galaxies {Dl, D2, D3} lose all 
of their halos and most of the stars in 6, 4.5, 
and 5 Gyr, respectively. As is the case with the 
large spirals, galaxy D2 from Cluster II experi- 
ences the strongest tidal heating, while galaxy Dl 
from Cluster I the weakest. The final state of the 
models is very peculiar. The particles are squeezed 
into a narrow "tube" stretching for several hun- 
dred kiloparsecs. The shape of the structure is 
truly one-dimensional; the axis ratios are smaller 
than a few percent. The orientation of the stream 
must be in the direction of the strongest compo- 
nent of the tidal force at the moment of disruption. 
Of course, such a long straight line is an artifact of 
the tidal approximation: in real clusters the tidal 
stream would bend and curve beyond a few tens 
of kpc. 

These results are not inconsistent with the pres- 
ence of dwarf galaxies in clusters. A continuous 
infall of small galaxies is expected. However once 
in the cluster, the diffuse dwarf sphcroidals should 
be disrupted in several crossing times. In contrast, 
the dense dwarf ellipticals should be almost unaf- 
fected by tidal heating and might be found in any 
parts of the cluster. 

6.2. Low Surface Brightness Galaxies 

A much more luminous but even more diffuse 
type is the low surface brightness galaxies. LSB 
galaxies have slowly rising rotation curves and 
very extended disks compared to the high surface 
brightness spirals. McGaugh & de Blok (1998) 
find i? ~ 10 kpc as a representative scale length of 
galaxies in their sample. Nevertheless, the asymp- 
totic values of the rotation curve of the high and 
low surface brightness galaxies are similar, and 
they both follow the same Tully-Fisher relation 
(Dalcanton et al. 1997b, a). 

Similarly to dwarf sphcroidals, the LSB galaxies 
are extremely sensitive to tidal heating. In order 
to investigate if they are also prone to destruction, 
I run three models {LSB1, LSB2, LSB3} varying 
the disk mass and scale length (Table 1). Each 
model is embedded in a massive halo. In model 
LSB1 the disk and the halo have the same masses 
as in galaxy L3 but spread out with i? = 10 kpc. 



In model LSB2 the disk mass is reduced by a factor 
of 4 relative to LSB1. And in the extreme model 
LSB3 the scale length is further increased to Rq = 
15 kpc. The tidal history of all LSB models is the 
same as in galaxy L3. 

Figure 12 shows the initial rotation curves for 
the three LSB models. Unlike Fig. 1, the x- 
axis scale is in kpc to emphasize the huge ex- 
tent of these galaxies. Even though the halo mass 
is the same in all three models, the maximum 
value of the circular velocity declines from 229 
km s _1 to 181 km s _1 as a result of the decreas- 
ing disk mass and increasing scale length. Dark 
matter dominates everywhere and determines the 
stability of stellar disks against tidal disruption. 
For galaxies LSB1 and LSB2 the average halo 
density in the core, 6 x 10~ 3 M Q pc -3 , is above 
the critical tidal density, but galaxy LSB3 with 
Pav(Rc) = 1-7 x 10~ 3 M Q pc~ 3 is not expected to 
survive. 

Simulations again confirm the predictions based 
on the tidal density. Figure 13 shows that the ha- 
los of models LSB1 and LSB2 are truncated at 
about 30 kpc, similarly to model L3. For the LSB 
galaxies, however, this implies losing 70% of their 
dark mass (see Table 3). As a result, the circu- 
lar velocity falls by 18% to 25%. Stellar disks 
arc transformed even more strongly. Although 
only 20% of the stars are lost, the rest settle into 
a spheroidal configuration with comparable axis 
ratios in all directions. The radial scale length 
shrinks, while the vertical scale height expands by 
more than a factor of five. The density distribu- 
tion is no longer exponential and declines more 
sharply (Figure 14). Any gas remaining in the 
diffuse disks would be stripped by ram pressure of 
the intracluster media. The likely outcome of such 
evolution is a gas-devoid spheroidal galaxy. 

The third model, LSB3, is essentially disrupted 
by tidal forces. Only 4% of the halo remains bound 
at the end of the simulation. Because of the larger 
initial extent and lower density, the truncation ra- 
dius is correspondingly smaller: R t ~ 30/1.5 = 20 
kpc (cf. Figure 13). The central density of the halo 
is reduced by more than a factor of four. Figure 14 
shows that the stellar component is also severely 
damaged. The disk shrinks to less than a half of 
its original size and loses most of the stars. The 
surface density has fallen at all radii. This small 
remnant of the galaxy is unlikely to survive tidal 
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heating for much longer and will be torn apart by 
future tidal shocks. 

To aid the comparison with observations, the 
critical tidal density (eq. [1]) can be translated to 
the projected surface brightness of the disk that 
would be disrupted by tidal forces. The face-on 
disks with the scale length Rq ~ 3 kpc would sat- 
isfy the disruption condition p av (Ro) < Ptid if their 
central surface brightness in the B band is below 

(31) 

or /i t id = 23.8 mag arcsec 2 in the R band. For 
the edge-on disks the critical surface brightness 
can be up to a factor of 10 higher. However, if the 
dark matter cores are even denser (ph,av > Pd,av) 
then nnd would correspondingly decrease. The 
discoveries in the Coma cluster of the giant dif- 
fuse arc with the surface brightness ps < 26.5 
mag arcsec -2 by Trentham & Mobasher (1998) 
and the diffuse plume-like feature with w 26 
mag arcsec -2 by Gregg & West (1998) may pro- 
vide the examples of a complete tidal disruption. 

6.3. Origin of Intracluster Light 

Diffuse emission, not associated with galaxies, 
contributes 10% to 50% of the light from clus- 
ters (Burstein et al. 1995; Vilchez-Gomcz 1999). 
As the tidal forces remove a substantial amount 
of stars from the low surface brightness galaxies, 
tidal streams stretch across the clusters. Is there 
enough LSB galaxies to provide the origin of dif- 
fuse intracluster light? 

The fraction of stars stripped from the galaxies 
can be estimated as follows. Assume that the ini- 
tial galaxy luminosity function is given by equa- 
tion (2) and that the luminosity, along with the 
mass, scales linearly with the initial extent of the 
galaxy: L oc R. Tidal truncation removes stars 
beyond the radius Rt ~ Rt,* (cg/cg,*) 5 where Rt,* 
is the tidal radius of L* galaxies. Then the lu- 
minosity after truncation is L t = L (R t / R) = 
L* (Rt.*/R*) {vg/vg,*), an d the fraction of dis- 
persed light as a function of x = L / L* is 

where the total luminosity of the cluster is L tot — 
J N(L) LdL = N L* T(2 - a). Using the adopted 



luminosity function and the Faber- Jackson rela- 
tion a g — <7 Sj * (L/L*) 1 / 4 , the distribution becomes 

The simulations show that i? t * w 30 kpc and the 
typical ratio R t ,*/R* = 0.2 — 0.5. For a = 1 — 1.5, 
the distribution peaks around the characteristic 
scale x = 1. 

This calculation shows that most of the stripped 
I s 1 ^s a c& £ m < e from giant galaxies. Since the gi- 
ants dominate the integrated light of both high 
and low surface brightness galaxies, the amount 
of the intracluster light in the tidal disruption 
scenario is proportional to the average fraction of 
LSB galaxies. The comparison of the APM survey 
and the CfA Redshift Survey (Sprayberry et al. 
1997) shows that LSB galaxies contribute about 
20% of the total light in the field. Dalcanton et al. 
(1997a) find that this fraction may even be higher. 
Also, in the Fornax cluster the LSB galaxies con- 
tain 20% to 30% of the total light (Kambas et al. 
2000). Thus the disruption of large low surface 
brightness galaxies would naturally account for 
the intracluster light. 

7. Comparison with Previous Work 

The idea of transforming spirals into SOs by 
tidal heating in clusters is not new. From the 
study of galaxy collisions in cluster cores, Rich- 
stone (1976) finds that the dark matter halos are 
severely truncated. Simulations of the tidal heat- 
ing of colliding galaxies (Farouki & Shapiro 1981) 
suggest that the disks grow in scale height, up 
to zo w 0.2 R . Despite the limitations of ana- 
lytical cluster models and small number of parti- 
cles, these early studies agree qualitatively with 
the more advanced computations. 

A detailed investigation of galaxies in clusters 
has not been possible until very recently. Moore 
et al. (1996) and Moore, Lake & Katz (1998) have 
conducted the first high-resolution simulation of 
the tidal effects, using a singular isothermal clus- 
ter model populated by a large number of tidally- 
truncated galaxies. Both the neighbor galaxies 
and the cluster halo exerted tides on the galaxy 
simulated with higher resolution. Moore et al. 
predicted a strong evolution: the disks of medium 
size spirals (V c — 160 km s -1 ) were destroyed 
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and transformed into dwarf spheroidals. Two ef- 
fects may have contributed to this: (i) they used 
the cluster with a higher velocity dispersion (like 
Coma) where the tidal forces are stronger than in 
the clusters (like Virgo) simulated here, and (ii) 
the clusters growing hierarchically had less time 
to exert tidal effects than a steady massive clus- 
ter. 

More recent simulations of Moore et al. (1999) 
have incorporated the hierarchical cluster forma- 
tion in the ilo = 1 cosmology and indicate more 
moderate tidal evolution. They find that the disks 
of high surface brightness spirals survive within 
their dark matter halos but expand in vertical di- 
rection by a factor 2 to 3. Their new results are in 
a much better agreement with our models L1-L3. 

Due to a smaller number of particles in the disk 
(2 x 10 4 ) and in the halo (5 x 10 4 ), Moore et al. 
(1999) simulations could suffer from numerical re- 
laxation more strongly than our results. Estimat- 
ing m h = 4 x 1O 6 M , Ph (R c ) w 0.2 Mq pc~ 3 , 
«cnc,h ~ 180 km s _1 for their isolated HSB galaxy, 
I find the disk heating time £ re i,h ~ 1-5 Gyr (cf. 
eq. [25]). However, their Fig. 5 shows that in 5 
Gyr the disk has expanded by only slightly more 
than a factor of 2, which corresponds to a relax- 
ation time of 4.5 Gyr. Thus for some reason, in 
Moore's simulation the numerical relaxation rate 
is three times slower. The only worrying sign is 
the steepening of the central density profile of the 
LSB galaxy in their Fig. 7. Since tidal heating 
vanishes in the center it cannot be responsible for 
this effect, which could instead be a manifestation 
of the artificial "core collapse" . 

Moore et al. (1999) have examined the survival 
of the disks as a function of their mass, scale 
length, and the core radius of the surrounding 
dark matter halo. They find that compact disks 
are stable regardless of the halo core radius and 
that extended disks are stabilized by denser ha- 
los. All this can be explained in terms of a single 
tidal density parameter (cf. eq. [1]). The galaxies 
are stable if their combined inner density, of stars 
and dark matter, is above the tidal density. The 
results of Table 1 of Moore et al. (1999) indicate 
the truncation radius of about 20 kpc for small 
halo core radii (R c = 1 kpc) and about 10 kpc for 
large core radii (R c = 10 kpc), all of which cor- 
respond to a tidal density p t id ~ 3.6 x 10~ 3 M 
pc~ 3 . (The fraction of stars lost is independent of 



the disk mass as long as the dark halo dominates, 
M d < M h .) 

Despite different cluster models and simulation 
methods, the critical tidal densities in Moore et 
al. (1999) and in this paper agree within a factor 
of two, and the predicted dynamical evolution of 
galaxies is also similar. Some of the difference may 
be due to the close encounters of neighbor galaxies 
within 60 h^ 1 kpc, which are not resolved in my 
cluster simulations (Paper I estimates that these 
encounters contribute 10% to 50% of the overall 
tidal heating, 7 t id, eq. [26]). Nevertheless, it is re- 
markable that such a simple concept as the tidal 
density of the Roche lobe provides a quantitatively 
accurate description of detailed simulations. It of- 
fers a direct way to compare the results of different 
authors. 

Tidal disruption of galaxies as the origin of in- 
tracluster light has been envisioned by Merritt 
(1984). He argued that the most probable time for 
disruption is during the cluster collapse. Moore 
et al. (1999) also confirm the disruption of LSB 
galaxies in their simulations. Calcaneo-Roldan et 
al. (2000) present an excellent example of a dif- 
fuse arc structure in the Centaurus cluster and 
explain it with the debris of a low surface bright- 
ness galaxy disrupted by the giant elliptical NGC 
4709. Their finding that the stripped stars pro- 
duce diffuse features like those observed goes along 
with the results of §6. On a scale of the Local 
Group, Mayer et al. (2001) find that the observed 
dwarf spheroidal galaxies could be the remnants 
of tidally-disrupted LSB progenitors. 

It should also be mentioned that Dubinski 
(1998) has combined the self-consistent cluster 
evolution with the galactic simulation in a single 
parallel tree-code run. He aimed to study the for- 
mation of a giant cD galaxy in the cluster center 
through a sequence of mergers of smaller galaxies. 
Unfortunately, his approach does not allow enough 
resolution to follow the dynamics of galaxies other 
than the central giant. 

8. Discussion 

The main result of this work is a natural and 
quantitative explanation of the transformation of 
spiral galaxies into SOs. The simulations show 
that even on the orbits with maximum expected 
tidal heating, large disk galaxies still remain flat- 
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tened and not elliptical. The evolution of elliptical 
galaxies must have proceeded prior to the cluster 
formation. 

The main advance is made with the supporting 
simulations of the hierarchical cluster formation, 
reported in Paper I. The galactic orbits, the tidal 
field along them, and the timcscalcs of dynamical 
evolution are all representative of the three plausi- 
ble cosmological scenarios of structure formation. 
The evolution is stronger in the models (£1 < 1) 
where galaxies spend more time in the dense clus- 
ter environment. 

Unlike previous studies with analytical cluster 
models, this work demonstrates that tidal heat- 
ing is not limited to the cluster core because large 
local substructures are present in all parts of the 
forming clusters. In fact, most of the observable 
effects (halo truncation and disk heating) happen 
in the outer regions, as soon as the galaxy enters 
the cluster. Also, most of the evolution can be pre- 
dicted using a single parameter, the critical tidal 
density along the orbit (eq. [1]). 

A testable prediction of the tidally-driven evo- 
lution is the preferential abundance of SO galax- 
ies in the middle of large clusters at the expense 
of the spirals. Disk galaxies on the cluster pe- 
riphery should have distorted profiles as well as 
the occasional spots of diminishing star formation. 
The abundance of elliptical galaxies, on the other 
hand, should be universal. These predictions are 
confirmed by observations (Jones, Smail & Couch 
2000). A systematic study of galaxies at various 
stages of their evolution in clusters by the Sloan 
Digital Sky Survey (Blanton ct al. 2002) will fur- 
ther improve the observational picture. 

New HST images of the cluster CL 1358+62 
(van Dokkum et al. 1998) show that the SOs are 
systematically bluer in the outer regions and may 
have formed stars more recently than the inner 
SOs. Interestingly, the brightest blue galaxies are 
bulge-dominated. Even though SOs are systemati- 
cally less luminous than both ellipticals and spirals 
(van den Bergh 1998), their bulges may appear 
more prominent after star formation has ceased: 
the disks would dim by up to 2 mag, while the 
bulges would brighten due to the reduced extinc- 
tion. These effects, however, may not be sufficient 
to explain the large bulges of SOs from the small 
bulges of Sa-Sc galaxies (Kent 1985). Therefore, 
like the elliptical galaxies, most of the bulges had 



to accumulate at the early stages of galaxy forma- 
tion. 

The most significant limitation of the present 
study is the absence of gas dynamics. Future work 
should include a self-consistent treatment of the 
gas, both in the cluster and the galaxy simulations. 
Parallel Tree-SPH codes should be able to address 
this issue, although the amount of computations 
is very large. Also, high-resolution cluster sim- 
ulations should help to resolve close galactic en- 
counters. A nested grid method is necessary here, 
since the dynamical range should be extended by 
at least a factor of 10, from 60 hr 1 kpc to less than 
5 h^ 1 kpc. Adaptive mesh codes, such as AMR 
(Bryan et al. 1995), could reach the desired level 
but the computational expense is again not cur- 
rently possible. The interpenetrating galaxy en- 
counters may be responsible for the unexpected 
star formation episodes, such as in the "E+A" 
galaxies (Gunn 1989; Zabludoff 1999). 

9. Conclusions 

I have used the high-resolution TV-body simu- 
lations to explore the dynamical evolution of disk 
galaxies in clusters. The evolution is driven by 
the tidal forces of the hierarchically forming clus- 
ters, extracted from the lower resolution simula- 
tions (Gnedin 2003) in three cosmological scenar- 
ios. I find that tidal effects are capable of induc- 
ing strong dynamical and morphological transfor- 
mations, which can be parametrized by the criti- 
cal tidal density along the galactic orbit. In the 
Virgo-type clusters this tidal density is p t id ~ 
2 x 10- 3 M Q pc- 3 . 

(i) Disks of large spiral galaxies thicken by a 
factor of two to three, while the stellar velocity 
dispersion in outer regions rises to 100 km s _1 . 
Toomre's Q parameter increases to > 2 and even 
a large amount of cold gas is unable to make the 
disks gravitationally unstable. The gas in the 
outer regions can be stripped by ram pressure of 
the intracluster gas. However, the stellar surface 
density distribution remains similar to the initial. 
Devoid of star formation these galaxies transform 
into the SOs, in agreement with the observations 
of low redshift clusters. 

(ii) Dark matter halos are truncated at 30 ± 6 
kpc, for the initial circular velocity of 250 km s _1 . 
In smaller halos the truncation radius is propor- 
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tionately smaller, R t cx V c . As a result, the galax- 
ies with extended halos may lose most of their 
mass to the cluster. 

(iii) Low density (dwarf spheroidal and low 
surface brightness) galaxies entering the cluster 
may get completely disrupted by tidal shocks. 
The stripped material would stretch in long tidal 
streams along the original trajectories. According 
to the current faint galaxy surveys, the unbound 
stars of giant LSB galaxies may provide enough 
material for the diffuse intraclustcr light. 

(iv) In all three cosmological models tidal heat- 
ing is not limited to the cluster core and starts to 
affect galaxies as soon as they enter the cluster. 
The variation of the evolution on different orbits 
within the same cluster is as large as the variation 
due to the different cluster models. Nevertheless, 
the epoch of cluster formation, the critical tidal 
densities, and the heating rates indicate that the 
low f2o models produce stronger dynamical evolu- 
tion of galaxies and agree with observations better 
than the flo = 1 model. 
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Fig. 1. — Initial rotation curve of large disk galax- 
ies (solid line) and the individual contributions 
from the exponential disk (dashes) and the isother- 
mal halo (dots). Cylindrical radius is in units of 
the disk scale length, i? - 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.0. 
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Fig. 2. — Upper panels: Initial conditions for large disk galaxies, face-on and edge-on views (only stellar 
particles are shown). Lower panels: Final configuration of galaxy LI. 
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Fig. 3. — Halo density profiles of large galaxies 
LI (dashes), L2 (dots), and L3 (solid line) at the 
end of each simulation. Outer thin line shows the 
initial distribution. 
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Fig. 4. — Rotation curves, with the disk and halo 
contributions, of large galaxies LI (dashes), L2 
(dots), and L3 (solid lines) at the end of each sim- 
ulation. 



Fig. 5. — Disk surface density of large galaxies LI 
(dashes), L2 (dots), and L3 (solid line) at the end 
of each simulation. Straight thin line shows the 
initial exponential profile. 




R/R 

Fig. 6. — Disk scale height, zo/Ro, as a function of 
cylindrical radius for large galaxies LI (dashes), L2 
(dots), and L3 (solid line) at the end of each simu- 
lation. Declining dotted line shows the numerical 
heating in the isolated run. Solid horizontal line 
is the initial scale height, 0.1 Rq. 
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Fig. 7. — Three components of the velocity dis- 
persion, a R, (70, o~ z , and the rotation speed V m t of 
the disks of large galaxies LI (dashes), L2 (dots), 
and L3 (solid line) at the end of each simulation. 
Dash-dotted line is for the isolated run. Thin solid 
lines in the three dispersion panels show the initial 
state of the disk. 
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Fig. 8. — The disk and halo masses of large galax- 
ies LI (dashes), L2 (dots), and L3 (solid line) 
versus time. Disk masses are essentially constant 
(horizontal lines). Halo masses drop steadily (de- 
clining lines), although the instantaneous amount 
of bound mass may fluctuate with the external po- 
tential. Note that simulation times differ for each 
cosmological model. Sudden mass losses can be 
identified with the peaks of the external tidal force 
(cf. Figs. 14-16 in Paper I): at t = 0.3,0.8,1.2 
Gyr for LI, at t = 1.8,2.2,2.8,6 Gyr for L2, and 
at t = 0.1, 1.4, 3.3, 4.2, 5.1 Gyr for L3. 
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Fig. 9. — Variation of the disk scale length with 
time for large galaxies LI, L2, L3 and for the iso- 
lated run. 
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Fig. 11. — Toomre's Q parameter for the disks 
of large galaxies LI (dashes), L2 (dots), L3 (solid 
line) at the end of each simulation. Lower solid 
line shows the initial state of the disk. 




Fig. 10. — Variation of the global scale height with 
time for large galaxies LI, L2, L3 and for the iso- 
lated run. In the three lower panels, the instanta- 
neous value of the scale height is normalized to the 
initial value and corrected for the effect of numer- 
ical heating (eq. [13]). Small dots show for com- 
parison the results of gadget simulations with 
more accurate vertical forces but with = 10 5 , 
N h = 4 x 10 5 particles. 
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Fig. 12. — Initial rotational curves, with the disk 
and halo contributions, of the low surface bright- 
ness galaxies LSB1 (dashes), LSB2 (dots), LSB3 
(solid line). 
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Fig. 13. — Final halo density profiles of galax- 
ies LSB1 (dashes), LSB2 (dots), LSB3 (thick solid 
line). Thin solid lines show the initial profiles for 
each galaxy. 




Fig. 14. — Final disk surface density of galaxies 
LSB1 (dashes) and LSB3 (thick solid line). Thin 
solid lines show the initial profiles for each galaxy. 
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Table 1 

Initial parameters of simulated galaxies 





M h 


M d 


V c 


Ro 




Trot 


Model 


(M ) 


(M & ) 


(km s -1 ) 


(kpc) 


(kpc) 


(yr) 


L1,L2,L3 


8 x 10 11 


4 x 10 10 


250 


3 


0.3 


1.2 x 10 8 


D1,D2,D3 


2 x 10 9 


10 8 


20 


1 


0.1 


5.0 x 10 8 


LSB1 


8 x 10 11 


4 x 10 10 


229 


10 


1 


4.6 x 10 8 


LSB2 


8 x 10 11 


10 10 


222 


10 


1 


4.7 x 10 8 


LSB3 


8 x 10 11 


10 10 


181 


15 


1.5 


7.0 x 10 8 



Table 2a 
Isolated runs 



Model 


N d 


N h 


At (yr) 


Ro/ Ro,i 




b/a 


c/a 


L-isol 


10 6 


10 6 


1.2 x 10 6 


1.04 


1.7 


0.99 


0.97 


L-isoll 


10 6 


10 6 


2.4 x 10 6 


1.06 


2.1 


1.00 


0.97 


L-isol2 


5 x 10 5 


5 x 10 5 


1.2 x 10 6 


1.04 


1.9 


1.00 


0.97 


L-isol3 


5 x 10 5 


10 6 


1.2 x 10 6 


1.04 


1.9 


0.99 


0.97 


LSB-isol3 


10 6 


10 6 


7.0 x 10 6 


1.01 


0.8 


1.00 


0.88 



Table 2b 
Energy changes in isolated runs 



Model 




(^E/Ef) d 


((AE/E) 2 ) h 


(A^) / (E z ) d 


L-isol 


1.3 x 10~ 3 


4.8 x 10~ 3 


1.3 x 10" 3 


1.0 


L-isoll 


-1.5 x 10~ 2 


1.1 x 10~ 2 


2.1 x 10" 3 


1.1 


L-isol2 


2.6 x 10~ 3 


7.2 x 10" 3 


1.7 x 10" 3 


1.1 


L-isol3 


4.9 x lO" 3 


7.9 x 10~ 3 


1.6 x 10~ 3 


1.2 


LSB-isol3 


3.2 x 10~ 2 


3.1 x 10" 4 


1.2 x 10~ 3 


0.55 
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Table 3 

Final parameters of large galaxies 



Model 


T sim (Gyr) 


M h /M h>i 


M d /M dti 


v c /v C}i 


Rol Ro,i 


z o/ z o,i 


b/a 


c/a 


LI 


9.35 


0.63 


1.000 


0.98 


1.01 


1.6 


0.93 


0.88 


L2 


10.71 


0.36 


0.995 


0.94 


1.04 


2.4 


0.93 


0.88 


L3 


12.28 


0.48 


0.999 


0.95 


1.03 


1.7 


0.89 


0.84 


L3a 


12.28 


0.55 


1.000 


0.99 


1.01 


1.4 


0.92 


0.87 


L3b 


12.28 


0.39 


0.997 


0.93 


1.02 


2.5 


0.91 


0.88 


LSB1 


12.28 


0.34 


0.81 


0.82 


0.74 


5.7 


0.86 


0.84 


LSB2 


12.28 


0.32 


0.77 


0.75 


0.78 


5.2 


0.87 


0.79 


LSB3 


12.28 


0.04 


0.28 


0.29 


0.43 


2.4 


0.79 


0.72 



Table 4 
Test models of galaxy L3 



Model N d 


N h 


At (yr) 


M h /M h ,i 


M d /M dti 


V c /V c ,i 


Ro/Ro,i 




b/a 


c/a 


main 10 6 


10 6 


1.2 x 10 6 


0.48 


0.999 


0.95 


1.03 


1.7 


0.89 


0.84 


1 10 6 


10 6 


2.4 x 10 6 


0.49 


0.998 


0.95 


1.08 


1.9 


0.89 


0.85 


2 5 x 10 5 


5 x 10 5 


1.2 x 10 6 


0.48 


0.999 


0.95 


1.04 


1.7 


0.89 


0.85 


3 5 x 10 5 


10 6 


1.2 x 10 6 


0.48 


0.999 


0.95 


1.04 


1.5 


0.89 


0.84 



Table 5 
gadget runs 



Model 


N d 




N h 


M h /M hti 


M d /M dti 


V c /V c ,i 


Ro/Ro,i 


zo/z ,i 


b/a 


c/a 


Ll-GADGET 


10 5 




10 5 


0.63 


0.999 


0.90 


1.15 


2.1 


0.90 


0.89 


L2-GADGET 


10 5 




10 5 


0.36 


0.994 


0.84 


1.10 


3.1 


0.94 


0.87 


L3-GADGET 


10 5 




10 5 


0.45 


0.999 


0.83 


1.18 


2.8 


0.93 


0.91 


Ll-GADGET 


10 b 


4 


x 10 5 


0.63 


0.999 


0.90 


1.12 


1.8 


0.90 


0.89 


L2-GADGET 


10 5 


4 


x 10 5 


0.36 


0.994 


0.84 


1.11 


2.5 


0.95 


0.88 


L3-GADGET 


10 5 


4 


x 10 5 


0.45 


0.999 


0.83 


1.20 


1.9 


0.92 


0.90 
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